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Abstract 

In the large iV c limit, the mass spectrum of the L = 1 orbitally excited baryons N* has a 
very simple structure, with states degenerate in pairs of spins J = (i, |), (|, | ), corresponding to 
irreducible representations (towers) of the contracted 5t7(4) c symmetry group. The mixing angles 
are completely determined in this limit. Using a mass operator approach, we study 1/N C corrections 
to this picture, pointing out a four-fold ambiguity in the correspondence of the observed baryons 
with the large N c states. For each of the four possible assignments, we fit the coefficients of the 
quark operators contributing to the mass spectrum to 0(N~ 1 ). We comment on the implications 
of our results for the constituent quark model description of these states. 
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I. INTRODUCTION 



Ever since their discovery, the negative-parity excited baryons have been a testing ground 
for the quark model, as their complex mass spectrum makes it possible to distinguish among 
different models of quark-quark interactions. In the quark model, the L = 1 negative parity 
baryons fall into the 70 representation of SU(6), which contains the (1,2), (10,2), (8,2) 
and (8,4) of SU(3) a ® SU(2) sp QQ. The finer details of the spectrum were first studied in 
a constituent quark model by Isgur and Karl 0,0]. Recently, the flavor dependence of the 
quark spin interactions has been a matter of some debate, related to their dynamical origin 
as arising from constituent gluons or pions exchange 0, |y, Q • 

In this paper we focus on a different, model-independent de scrip tion of these states, 
which was proposed only relative ly r ecently 0, EH EH EH EH EH Il7l EH EH HH Hl| . This 
is based on the large N r limit [H, HI, and makes use of the contracted SU(2Nf) symmetry 
of QCD in this limit There are two equivalent ways of extracting the predictions of 

this symmetry. The first method is algebraic, and makes use of commutation relations of 
operators representing physical quantities such as axial currents and hadron masses, which 
are constrained by so-called consistency conditions. Although somewhat cumbersome and 
difficult to extend beyond the leading order in l/N c , such an approach has been applied to 
the study of the mass spectrum and strong decays of the excited light baryons in [IH Efif . 
These states were predicted to fall into irreducible representations of the contracted SU(2Nf) 
algebra, with nontrivial implications for the mass spectrum and mixing angles of the excited 
baryons. 

A somewhat different approach to large N c baryons is based on the quark operators 
method. In this approach the representations of the contracted SU(2Nf) algebra are con- 
structed in terms of a 'quark' basis ||, S3) (g) If,^). Although similar to the one-quark basis 
states used in the constituent quark model (CQM), we stress that this is a purely mathemat- 
ical device and does not imply any of the dynamical assumptions of the CQM. Any physical 
quantity is further written as an expansion in n-body operators acting on the 'quark' states, 
ordered according to n. This can be organized as an expansion in 1/N C , with a finite number 
of independent operators contributing at any given order in 1/N C . This method for extract- 
ing the predictions of the large N c limit has the advantage of being more systematic and 
easier to extend to subleading orders in 1/N C . 

The mass spectrum of the excited baryons has been studied using this approach in a 
series of papers HI El IH |2l|, including operators of order up to and including 1/N%. In 
these papers the coefficients of these operators have been extracted from a fit to the observed 
masses, yielding a very specific pattern of the dominant operators. 

In this paper we reanalyze the 1/N C expansion for the masses of the nonstrange excited 
baryons. Guided by the expected structure of the large N c mass spectrum, we present fits 
for the coefficients of the mass operator, order by order in N c . We point out a four-fold 
ambiguity in the correspondence of the physical states with the large N c states, and the 
special status of the excited A* states in the N c expansion. At leading order in N c only 
three operators contribute to the mass operator, and we discuss the four equivalent fits to 
their coefficients, one of which turns out to be ruled out. Including l/N c corrections allows 
eight operators to contribute to the mass matrix. Considering only data on N* masses, the 
coefficients of these eight operators can be determined as functions of the two mixing angles, 
up to an additional continuous ambiguity. We present restrictive constraints on the allowed 
set of coefficients, first from dimensional analysis and then by including data on the excited 
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A* states. We compare the results of our fits to those previously discussed in the literature. 



II. FORMALISM 



The spin-flavor structure of the orbitally excited L = 1 baryons for arbitrary A c can be 
conveniently described in terms of quark quantum numbers. The quark basis states can be 
labeled as | J, /, S; mj, I3) where 5* is the quarks' spin, / their isospin and J = S + L the total 
baryon spin. The spin-isospin wave function of a nonstrange L = 1 orbitally excited baryon 
has mixed spin-flavor symmetry. This implies that for a given isospin /, the spin 5* takes all 
values compatible with \S — I\ < 1 except for S = I = N c /2. The most general state \JIS) 
can be constructed by adding one excited quark to a symmetric 'core' of quantum numbers 
\S C = I c ), such that S = S c + s, I = I c + i with s — i — |. 

In the large N c limit the mass eigenstates obtained from this construction fall into irre- 
ducible representations of the contracted SU(2Nf) algebra |1Q|. These are infinite towers 
of states, labelled by an integer or half-integer T, and contain all possible states satisfying 
the condition | J — I\ < T. All the states belonging to a tower are degenerate in mass, and 
the separation between towers is of order 0(N®) in the large A c limit. On the other hand, 
mass splittings within towers are small, and appear first at order 0(1/N C ). The tower states 
\JIT) are related to the quark model states \ JIS) by a simple recoupling relation [l5[ |24| 



| J/T) = ( _ ) m + L + ^ v / (2S + 1)(2T + 1) |^ y T ^\JIS). (1) 

This completely determines the mixing matrix of the L = 1 excited baryons in the large A c 
limit, in the absence of configuration mixing [Hf. 

The resulting mass spectrum for the nonstrange P— wave baryons is shown in Fig. 1(b). 
There are five 1 = 1/2 states (A— type) which correspond to quark model states of spin 
S = 1/2,3/2, and eight I = 3/2 states (A— type, not shown), corresponding to quark spin 
S = 1/2, 3/2 and 5/2. This covers all the observed A type states, but due to the finite value 
of A c = 3 in the real world, only two 1 = 3/2 states are present. This will be seen to limit 
the predictive power of the large A c approach for the A— type excited baryons. 

In the following we show how the large A c predictions are reproduced in the quark 
operators language. This will allow us to include also symmetry breaking to this picture, 
appearing at order 1/A C . The mass matrix of the L = 1 baryons can be written as a sum of 
operators acting on the quark basis as 

M = E j^C h O h (2) 

fc=0 JV c 

with Ok a fc-body operator. Both the coefficients and the matrix elements of the oper- 
ators on baryon states (Ok) have power expansions in 1/A C with coefficients determined by 
nonperturbative dynamics 

^ = E^t\ (o*) = E^W (n) - (3) 

n=0 Jv c n=0 Jv c 

The natural size for the coefficients C { k n) is A/A™, with A ~ 500 MeV. The basic building 
blocks for constructing the operators Ok are s\ r a , g ia (acting on the excited quark), Si, 
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(3/2, 5/2) 

N 3/2 (1700) T = 2 

N 5/2 (1675) 

N 1/2 (1650) 021 T = 



N 1/2 (1535) (1/2, 3/2) 

N 3/2 (1520) ' * 



(a) (b) 

FIG. 1: The mass spectrum of the nonstrange L = 1 orbitally excited nucleons. (a) the observed 
states and (b) in the large N c limit these states fall into three towers labelled by the quantum 
number T = 0, 1, 2. 

(acting on the core), and operators acting on the orbital degrees of freedom l 1 ,^ 2 ' 1 ^ = 
\{l\P} - \l 2 8» . 

We will use the operator basis introduced in which contains the following operators, 
ordered according to their leading power in 1/N C . At leading order in N c only three operators 
contribute to the mass matrix, given by 

O x = N c l , 2 = ZV , O s = J-Z^V^f . (4) 



One small technical difference with the operators in [18( is that we define Os with an added 
factor of 3. With this redefinition, its matrix elements are not anomalously small, which 
is required for a consistent dimensional analysis for its coefficients |2p| . l2f| . At subleading 
order (^(iV" 1 ) five additional operators start contributing, which can be chosen as 

Ot = ls + j^-^ltG e , 5 = ^lS c , O e = r^S c S c , (5) 

e> 7 = -Us c , o 8 = ^-i (2) ss c . 

N c N c 

These operators have a direct physical interpretation in the quark model in terms of one- 
and two-body quark-quark couplings. 

We start by examining the mass spectrum of the iV* states in the large N c limit. Keeping 
only the Oi,2,3 operators contributing at 0(N®), one finds by direct diagonalization of the 
mass matrix the mass eigenstates in the large N c limit as linear combinations of the quark 
model N\/2, N[, 2 states (apart from the redefinition of O3, we use everywhere the notations 
of 0) 

|r = o,j=i> = + (6) 

\T=l,J=\) = -f-N 1/2 + ±=Ni l/2 (7) 
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with masses 



M (0) = N c cf ] - Cf - -Cf ] (8) 

8 

= N C C[ 0) ~ icf + Aci°> . (9) 
A similar diagonalization of the mass matrix for the J = | iV* states gives the eigenstates 

|T=l,J=jb = ^3/ 2 + J^iV' (10) 



2 ; " 76 iV3/2 ^V6 ^ 

3 \ /5 Ar 1 

2> = -V6 JV ^ + ^ 



|r = 2,j = -> = -,/-jv 3/2 + -=jv^ 2 (ii) 



with masses 



Mf > = - Icf + Acf (12) 

Mf = tf e cf + icf -lef . (13) 

The A^ 5 /2 state does not mix and has the mass M^ "*. 

These results make the tower structure shown in Fig. 1(b) explicit, with four of the 
states degenerate in pairs (I, J) = (|, |), (|, |) (corresponding to the T = 1 tower) and 
(I, J) = (|,|))(|)|) ( corres P on ding to the T = 2 tower), respectively. Also, the mixing 
matrices agree with the expected recoupling relation ((TJ) and do not depend on the explicit 
values of the coefficients C^l, the reason being that O2 , O3 commute in the large iV c limit. 

The mass splittings between the towers are determined by the coefficients and C^°\ 

There is a discrete ambiguity in the correspondence of the large N c tower states and the 
observed iV* excited nucleons. Taking the number of colors N c to be infinitely large, the 
states arrange themselves into towers as described above. The observed mass values (see 
Fig. 1) suggest assigning the two lowest lying states iV(1520) and iV(1535) into the T = 1 
tower, and then group the states of spins 3/2 and 5/2 into the T = 2 tower. The remaining 
J = 1/2 state iV(1650) would belong to the T = tower. This was the assignment suggested 
in [3] as it comes closest to the large N c limit picture. 



Assignment 


T = 


T = 1 


T = 2 


# 1 
#2 
#3 
#4 


iV 1/2 (1650) 
iV 1/2 (1535) 
iV 1/2 (1535) 
A 1/2 (1650) 


{N 1/2 (1535), N 3/2 (1520)} 
{N 1/2 (1650), iV 3 /2 (1520)} 
{iV 1/2 (1650),iV3/2(1700)} 
{iV 1/2 (1535),iV3/2(1700)} 


{iV 3 /2(1700),iV5/2(1675)} 
{iV 3 / 2 (1700),iV5/2(1675)} 
{iV 3 /2(1520),iV5/2(1675)} 
{iV 3 / 2 (1520),iV5/2(1675)} 



Table 1. The four possible assignments of the observed nonstrange excited 
baryons into large N c towers with T = 0, 1, 2. 



However, this is not the most general possible assignment. Allowing for potentially large 
1/N C corrections, there are four possible assignments of the observed states into towers, 
defined in Table 1. 
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Each of these assignments leads to a different picture in the large N c limit, and to different 
predictions for the properties of these states. For example, the large N c predictions for ratios 
of strong decay amplitudes depend on the tower assignment of the states |l5j . In the next 
Section we determine the coefficients C[°\ of the operators appearing in (J2J) by imposing the 
constraint that in the large N c limit the physical states go over into the towers corresponding 
to each assignment. 

The status of the excited A* states appears to be special. In the large N c limit, there are 
eight 7 = 3/2 (A— type) states, which arise in the quark model from states with quark spin 
S = 1/2, 3/2 and 5/2. Coupling the spin S with the orbital momentum L = 1 gives 2 A^ 2 , 

3 Ag/ 2 , 2 Ag/ 2 and one Ay 2 states, where the subscript denotes the total J = S + L hadron 
spin. Diagonalizing the mass matrices of these states, one expects the mass eigenstates to 
arrange themselves again into towers, degenerate with the N* states, as follows: one Ag, 2 

state with mass corresponding to the T = tower, 3 states A*^ 2 3 / 2 5 , 2 degenerate with 

mass (the T = 1 tower), and 4 states A*/ 2 3 , 2 5 / 2 7 , 2 degenerate with mass M 2 (the 
T = 2 tower). 

The real world is very different: at N c = 3 only two A* states exist, with spins J = 
1/2, 3/2. Consider for example the mass matrix of the two A*^ 2 states for arbitrary N c > 5. 

Expressed in the basis of the quark model states | J IS) = ||§|) and ||||) it has the form 

/ M^NJ m 12 (N c )\ 
~ { m 12 (N c ) M 2 (N e ) ) ■ (14) 

The matrix elements Mi^ 2 (N c ) , m,i 2 (N c ) have power expansions in 1/N t 

M X {N C ) = N c d{N c ) + ~C 2 (iV c ) + 0{N- 1 ) , (15) 
M 2 (N C ) = N cCl (N c ) + 1J ^^C 2 (N C ) + 2N ' + ™ C ~ 15 C 3 (N C ) + O(N^) , (16) 



^iN c) = ^ (c 2i N c) - 3 -^C. iiNc) ) + O(iV-) . (17) 

Denoting the eigenvalues of this matrix with Xi t2 {N c ), they can be also written as expansions 
in 1/N C , X h2 (N c ) = N c d(N c ) + 0(N°) H . It is easy to check that the eigenvalues in 

the large N c limit are identical with the tower masses M{ 2 given in Q, (|T3*|) . At N c = 3 
the 'phantom' state S = 3/2 disappears, and its mixing with the 'physical' S — 1/2 state 
vanishes. The mass of the observed state A^ 2 is identified with the diagonal matrix element 
A/i(Ao — 3). 

The meaning of the large N c expansion for this case is less clear than for the iV* states, 
since the N c = 3 and large N c sets of states are completely different. In particular, it is 
not completely clear that the large N c expansion can predict the masses of the A* states 
with corrections parametrically suppressed by 1/N C , as for the iV* states. For example, the 
0(N°) term in the expansion of the A*^ 2 mass ()15J) differs from both tower masses Q, ()13|) 
(to which it goes as N c — > oo) by terms of 0(N®). 

For these reasons we will perform our analysis below in two steps. First, we exclude the 
A* states and explore the form of the most general constraints which can be obtained from 
the iV* alone. Then we add them in a second step, which allows us to compare our results 
with previous fits in the literature, where they are taken into account. 
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III. NUMERICAL RESULTS AND FITS 



The general form of the mass matrix with N c = 3 is rather complicated and involves 
mixing among states with the same quantum numbers. Therefore a direct analysis is not 
very transparent, and we present first a simplified fit at leading order in 1/N C . For this 
purpose, we keep only the 0(N c ,N®) terms in the mass matrix, which come from the unit 
operator 0\ and the leading terms in the expansion of the matrix elements of 02,3 • in a 
second step we allow for the contributions of all l/N c operators, which can reproduce the 
finer structure within the towers. 

We start by determining the values of the coefficients c[°l 3 in the large N c limit. The 
large N c mass eigenvalues and eigenstates have been presented in Eqs. © and (fT^j). 
For each assignment, we fitted the coefficients c[°l 3 to the observed N* masses. We use the 

experimental values for the masses given in Table V of This gives the results for C^2 3 
shown in Table 2 for each of the four possible assignments, together with the mixing angles 
@ni,@N3, which are fixed to the large N c values as given in ©, ©, (fTU)) . (JTTJ). We define 
the mixing angles as in |l3[ 



and 



iV(1535) \ 
A(1650) J 



iV(1520) \ 
A(1700) J 



cos 6*7vi sin#/vi 
— sin#/vi cos#ati 



COS sin#AT3 
-sin 6^3 cos 0^3 




N 3/2 

K,2 



(18) 



(19) 



As mentioned in Sec. II, the best fit is obtained for the assignment #1, for which the observed 
N* mass spectrum comes closer to the tower structure in Fig. 1(b). 

These results must satisfy an additional constraint, following from the no-crossing prop- 
erty of the eigenstates with the same quantum numbers. Consider for example the masses 
of the J = 1/2 states as functions of 1/N C . They can not cross when N c is taken from 3 to 
infinity, which means that the correspondence of the physical iVi/2 states with the large N c 
towers is fixed by the relative ordering of the T = 0, 1 towers. A similar argument can be 
given for the J = 3/2 states which fall into the T = 1, 2 towers. 



Assignment 


A c Cf 


cf 


CP 


X 2 /d.o.f. 


M 


M x 


M 2 


9m 


® N3 


#1 


1625 ±8 


83 ± 14 


-188 ± 28 


0.38 


1660 


1525 


1679 


2.53 


1.15 


#2 


1617±8 


115 ±14 


-57 ±27 


20.3 


1538 


1542 


1679 


0.96 


1.15 


#3 


1615 ±9 


-12 ±16 


142 ± 38 


94.1 


1538 


1666 


1601 


0.96 


2.72 


#4 


1592 ±9 


2± 16 


-111 ±40 


98.5 


1660 


1557 


1601 


2.53 


2.72 



Table 2. The values of the coefficients C[°} 3 (in MeV), the masses of the towers 
^o,i-2 (in MeV) and the mixing angles (in rad) following from the large N c fit, 
for each of the 4 assignments. 

Taken together with the ordering of the observed hadron masses, this gives the following 
correspondence of the four assignments with the ordering of the tower masses: 



#1 



{M ,M 2 } > M 



(20) 
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C2 



FIG. 2: Regions in the (C 2 , Cg ) plane corresponding to each assignment. Each area corresponds 
to one assignment: #1 (black), #2 (dark grey), #3 (light grey), #4 (white). The line separating 
the assignments #3, #4 and #2,#1 is given by 3C3 + gC^ = 0, and the line separating the 
regions for assignments #2,^3 and #1,#4 is 3C^ — 8C^ = 0. The strips show the results for 
the coefficients C23 following from the large N c fit in Table 2, for each assignment. 



#2 : M 2 >M 1 > M , 
#3: M 1 >{M ,M 2 }, 
#4 : M > Mi > M 2 . 

Furthermore, since the relative ordering of the tower masses is given by C^l, there is a unique 
assignment corresponding to each set of these coefficients. This correspondence is shown in 
a graphical form in Fig. 2, together with the results for the coefficients C 2 °l following from 
the fits to the physical masses in Table 2. 

The plot in Fig. 2 can be used to further constrain the values of the coefficients C23 ■ Note 
that the region for these coefficients obtained for assignment #2 crosses partially into the 
region corresponding to This restricts further the range of values for C^l corresponding 
to #2. Also, the fit results corresponding to #4 has no overlap with the area corresponding 
to this assignment, which means that the assignment #4 is ruled out at leading order in N c . 
The results of this analysis can be compared directly with a 3-parameter fit performed in 
. where only the operators ©1,2,3 were included (see Table V in (HI)- The only differences 
are the inclusion of the A* states in 0] and the fact that the authors of 0] use in this fit 
the matrix elements of the operators ©1,2,3 at N c = 3, while we evaluate them in the large 
N c limit. The coefficients obtained in are N C C X = 1626 ±6, C 2 = 75 ±9, C 3 = -146 ±17 
(in MeV), and come closest to our assignment #1. However, when all operators up to order 
0(1/ N c ) are included, the coefficients obtained in (reproduced here in Table 3) come 
closest to our assignment #3. 
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N C C, C 2 


c 3 


c 4 


c 5 


c 6 


C 7 


C 8 


1389 ±60 -36 ±41 


123 ± 69 


87 ±97 


86±<< 


30 438 ± 102 


-40 ± 74 45 


± 172 



Table 3. Results for the coefficients Ck (in MeV) obtained from the fit of 18], 
using as input parameters the excited N, A* and mixing angles (9m , 9nz) = 
(0.61,3.04) (in rad). Note that our C 3 is related to that used in [18[ by C 3 = 
C 3 CCGL /3. 

Beyond leading order in 1/N C there are three types of symmetry breaking corrections 
which introduce deviations from the tower picture: 

1. Subleading corrections C[ 1 } 3 in the coefficients of the operators O1-3. 

2. 0(1/N C ) subleading corrections to the matrix elements of the operators 02-3- 

3. Leading contributions from the matrix elements of the O^s operators. 

While type 1. corrections just shift the relative position of the towers, the remaining two 
types of corrections introduce splittings among the tower states. We will include subleading 
corrections to the matrix elements of the operators Ois by evaluating them at the physical 
value N c = 3. 

Using as inputs the five masses of the N* states, it would appear at the first sight that 
there is a 3-dimensional manifold of solutions for the coefficients Ci__§ of the full set of 
0(1/N C ) operators. In fact this is reduced to a 2-dimensional set of solutions because of an 
accidental relation among the matrix elements of the mass operators with N c = 3. Namely, 
when restricted to the subspace of the N* states, a certain linear combination of Oqj turns 
out to be equivalent to the unit operator 

(0 6 ) No=3 -(0 7 ) Nc=3 =~-l. (21) 

This relation implies an ambiguity in the coefficients Ci^j beyond leading order in 1/N C , 
which can be summarized by noting that the mass spectrum of the N* states is invariant 
under the simultaneous changes 

N c d -> N c d + 5, C 7 ^C 7 + 25. (22) 

In other words, only the two combinations N c Cf s = N c Ci — \C 7 and C| ff = G% + C 7 can 
be determined unambiguously from the N* masses. Of course, the relation (|21j) is not true 
anymore when considering also the A* states, which can be used to resolve the ambiguity 
(|2*2*|) . However, as discussed in Sec. II, due to the incomplete structure of the towers in the 
1 = 3/2 sector at iV c = 3, the mass spectrum of the observed A* states does not have the 
correct large N c limiting behaviour. Therefore, we will keep our discussion as general as 
possible, and present separate results with and without including the A*s. 

To present the results of our numerical analysis for the coefficients Ci_g, we need to pick 
a parameterization of the 2-dimensional manifold of solutions. A particularly convenient 
choice of the coordinates on this manifold are the mixing angles (Onions), which will be 
left completely arbitrary. Using as inputs the masses of the five iV* states one finds, at each 
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0.5 1 1.5 2 2.5 3 

Si 



FIG. 3: The different regions in the (9ni,Gn3) plane (in rad) corresponding to the four assignments 
in Table I, determined as explained in the text. The black region - jf= 1, dark grey - # 2, light grey 
- # 3, white - # 4. The four dots shown correspond to the values of the mixing angles in the large 
iV c limit. 

point in the (9m, 9m) plane, one set of coefficients Cf|,C 2 _5, 8 . The analytical expressions 
for these coefficients can be found in the Appendix. 

It is natural to ask if there is a unique association of a particular solution at a given point 
(9m, 9ns) with one of the four assignments in Table 1. As noted above, the assignment which 
is realized depends on the large N c values of the cofficients C^l - see Fig. 2. Requiring that 
the extracted values of 6*2,3(^1,^3) lie inside the regions in Fig. 2 associated with each 
assignment, one finds the corresponding partition of the entire (9m, 9ns) plane into regions 
shown in Fig. 3. This procedure assumes that the 1/N C corrections to the coefficients Cfl 
are negligible; including them will change slightly the boundaries of the regions in Fig. 3. 

Without additional input, the individual coefficients C\s can not be determined from 
this 1/N C analysis, apart from the rather loose absolute bounds (obtained by varying 9m,N3 
over their entire range [0,tt)): 1374 < N c Cf < 1770, -148 < C 2 < 164, -326 < C 3 < 363, 
-169 < C 4 < 223, -169 < C 5 < 342, -251 < Cf < 463, -638 < C 8 < 775 (in MeV). 

A more precise determination becomes possible if we impose constraints on the parameter 
space from dimensional analysis. As discussed in Sec. II, the natural size of the coefficients 
is ~ A/iV™. In particular, this means that the coefficients 6*2,3 can not be too different 
from their large N c values extracted in the previous step (see Table 2). By simple power 
counting, one expects the 1/N C corrections to these coefficients to be of order A/N c ~ 150 
MeV. Restricting the coefficients C*i_ 3 singles out a solution in the (9m, 6W3) plane, resolving 
the ambiguity ([22)1 and fixing the remaining coefficients C*4_g. 

We present in Fig. 4 the regions allowed from such constraints, for each of the four 
assignments introduced in Table 1. The regions in the (9m, 9ns) plane are obtained from 
imposing the condition that the coefficients 6*2,3(^1,^3) do not differ by more than 150 
MeV from their large N c central values in Table 2. The results for the assignment #4 are 
given only for illustrative purposes, given that it is ruled out at leading order in N c . For 
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0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

9i Si 



(c) (d) 

FIG. 4: Regions in the (0jvi,0jv3) plane corresponding to (C2,Cs) close to the large iV c values in 
Table 2, for each assignment: (a) - # 1, (b) - # 2, (c) - # 3, (d) - # 4. The contour plots are given 
by max|C 2 , 3 - c£§| = 50, 100, 150 MeV. 

example, the black areas in Fig. 4(a) denote the region in the (9ni,@n3) plane for which 
both C*2 and C3 are within 50 MeV from their central values in Table 2, corresponding to the 
assignment #1. In the 6*2,3 plane (Fig. 2) this region is a square of sides 100 MeV centered 
on (83 MeV, —188 MeV). Each lighter shade of grey corresponds to an additional 50 MeV. 
Similar contour plots are given in Figs. 4(b)-4(d) for each of the remaining three assignments 
in Table 1. Although the allowed regions appear disjointed, they really form one continuous 
area because of the doubly periodic conditions in the (#jvi,#jV3) plane. 

It is not very illuminating to quote ranges of values for the coefficients Cf. correspond- 
ing to the allowed regions in Figs. 4(a)-(d), since there are significant correlations among 
coefficients. Instead we give in Table 4 their values at a few representative points in these 
plots, for each assignment. These points correspond to minimal 1/N C corrections to the 
coefficients C%1 as determined in Table 2 at leading order in N c . We quoted the coefficients 
Ci 5 6,7 in a form invariant under the ambiguity (J2*2*|) . One can choose to vary N C C\ about its 
large N c values in Table 2 over a range ~ A, which introduces a large uncertainty in the 
individual values of Ci,6,7- On the other hand, the combinations N C C[ S = N C C\ — \Cj and 
Cg = Cq + C-j have smaller errors, of order A/N r . 

Finally, to compare with the results of we consider also the constraints from the 
A* states. Although these states can not be included in a meaningful way in the large 
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FIG. 5: The constraints in the (Oni,0N3) plane following from the A*^ 2 and Ajy 2 masses. The 
black (grey) regions correspond to the Ay 2 — A*^ 2 = 75 ± 55 (75 ± 80) MeV. 



N c fit discussed at the beginning of this section, they can be used to constrain the 1/N C 
determination. Their masses are given explicitly by 



A* 



A 



3/2 



N C C X + hj 2 - jjc 5 + jjc 6 - l -C 7 



(23) 
(24) 



Their splitting Ag /2 - A*, 2 = -|C 2 + |C 5 = 75 ± 80 MeV gives an allowed region in 
the (0jvi,0jV3) plane, shown in Fig. 5. This constraint includes the point (0 n1 ,8 N3 )q M = 
(0.61, 3.04) obtained from the quark model fit in ^3|. While most of the region corresponding 
to the assignment #1 is excluded, this assignment is still marginally allowed (see point 1' in 
Fig. 5). Note that each point in the dark area of Fig. 5 gives an acceptable 0(1/N C ) mass 
matrix of the observed N*, A* states, and generalizes the fit of which corresponds to 
the point QM and is obtained by requiring agreement with the quark model mixing angles. 





0N1 QnZ 


N c Cf C 2 C 3 




#la 
#16 


2.18 1.18 
2.93 0.50 


1694 83 -187 
1446 84 -188 


-5 28 -114 -120 
-14 -97 333 -272 


#2a 
#26 


0.19 0.58 
2.72 1.34 


1466 116 -58 
1670 115 -58 


-29 -138 297 -46 
-125 -56 -71 205 


#3a 
#36 
QM 


1.51 1.77 
2.86 2.06 
0.61 3.04 


1758 -11 143 
1611 -13 142 
1410 -35 123 


41 211 -230 360 
-111 109 35 451 
85 84 398 49 


#4a 
#46 


1.32 0.37 
3.00 0.10 


1508 2 -110 
1379 2 -110 


186 97 222 -360 
22 454 -245 



Table 4. Coefficients (in MeV) and angles (in rad) for the points shown in 
Fig. 4(a)-(d). We also show the quark model point (QM) corresponding to the 
angles 9m,m used as inputs in j[18j . 
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IV. CONCLUSIONS 



We studied in this paper the mass spectrum of the L = 1 orbitally excited nonstrange 
baryons in the 1/N C expansion. In the quark model with spin-flavor independent quark 
forces, these states fall into the 70 representation of SU(6) [1]. While large N c QCD ar- 
guments can be used to reproduce many of the predictions of the quark model for these 
states in a model-independent way [15j| . they also predict a nontrivial structure of the mass 
spectrum which breaks spin-flavor symmetry [14j . Still, a significant amount of symmetry 
remains, in the form of the contracted SU(2Nf) symmetry of large N c QCD [1, This 
symmetry predicts that in the large N c limit, the L = 1 orbitally excited baryons arrange 
themselves into 3 infinite towers of degenerate states 2> 16 1. These towers are labeled by 
one integer T = 0, 1, 2, and contain all states which satisfy | J — I\ <T. 

In a parallel development, the structure of the mass spectrum of the excited baryons was 
studied in the 1/N C expansion 3 12 , H, 2^, 21 1 using the different but equivalent approach 
of quark operators. In this approach, the mass operator is written as the most general sum 
of flavor singlet Lorentz scalar operators which contribute to a given order in 1/N C . Using 
such method, the mass spectrum of the L — 1 baryons was studied in both 577(2) and 577(3) 
flavor, up to subleading order 0(1/ N^). 

In this paper we show explicitly how the predictions of the SU(4) C symmetry for non- 
strange excited baryons follow from the quark mass operator approach in the large N c limit. 
Apart from recovering the known structure of the mass spectrum in the symmetric limit, 
this approach allows also the study of the 1/N C symmetry breaking corrections. We discuss 
the special status of the excited A* states in the large N c expansion; for N c > 3 there are 8 
A* states (vs. only 2 at N c = 3), and their mixing has to be taken into account in order to 
reproduce the large N c symmetry predictions. This casts some doubt on the ability of the 
1/N C expansion to correctly describe their properties. 

The main aim of our paper is to answer the two (related) questions: a) are the predictions 
of large N c SU(2Nf) c symmetry still visible in the observed mass spectrum of the excited 
baryons? and b) to determine the values of the coefficients of the various operators in the 
mass operator. 

At the first sight, the structure of the observed mass spectrum of the N* states 
is very similar to the expected set of towers, with states almost degenerate in pairs 
{A^i /2 (1535), A" 3/2 (1520)}, and {A" 3/2 (1700), A" 5/2 (1675)}. This corresponds to a special as- 
signment of the physical states into SU(4) C representations (#1), and is supported by the 
full 1/N C analysis. This assignment allows for small coefficients for the 0(1/N C ) operators 
C*4_8 in the mass matrix (see Table 4, #la). Even after including the excited A* (see Fig. 5), 
it is marginally still allowed, in a region around (9m, 9m)i' ~ (2.0, 1.4). We conclude that 
the assignment #1 is still allowed within the present experimental uncertainties. 

Allowing for larger 1/N C corrections, 3 other assignments of the observed states into 
large N c towers are possible (see Table 1). We show that model- independent constraints in 
the large N c limit rule out one of them (#4). The remaining two assignments require (as 
expected) relatively large coefficients of the 0(N~ 1 ) operators (see Table 4), depending on 
the precise position in the (9m, 9n3) parameter space. Beyond leading order in N c there is a 
double continuum of coefficients Ck(9m, 9m) (with one additional ambiguity when excluding 
the A* states). We use dimensional analysis to study the regions in the parameter space 
where the 1/N C expansion is well behaved (see Figs. 4). 

We find that, despite requiring large 0(N~ 1 ) coefficients, the assignment #3 is currently 
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favored by the constraints from the A* mass splittings. This agrees with the results of the 
fit in [la ], where a special solution for the coefficients of the mass operator was determined 
using as inputs the masses and mixing angles of the orbitally excited N*, A* states. The 
specific pattern of the operator coefficients in this fit was invoked in the literature 0, 0| as 
supporting a model of interquark forces based on Goldstone boson exchange (as opposed 
to one-gluon exchange). Although clearly interesting, we find that it is premature to draw 
any such conclusions based on the large N c expansion. As noted above, even in the large 
N c limit, the mass operator can be determined only with discrete ambiguities; at subleading 
order this becomes a double (triple) continuous ambiguity, depending on whether the A* 
states are (are not) included. 

Further progress towards resolving this multiple ambiguity can be made by considering 
also the strong couplings of these states. In the large N c limit, these were shown to be 
related to the assignment of the physical states into large N c towers [l5| . Unfortunately, the 
precision of the present data is still not sufficient to resolve the ambiguity . Progress can 
also come from lattice QCD, from which first results on excited nucleons are already available 
22l | . As discussed in Sec. Ill, the specific assignment realized in Nature depends on the values 
of the large N c coefficients C23, which could be computed in a lattice simulation. Finally, 
improvements in the experimental uncertainties in the masses of these states (which will be 
available from Jefferson Lab in the next future) could sharpen the constraints presented here. 
Clearly, more than forty years after their discovery, the orbitally excited baryons remain an 
active and exciting field of study. 

Note added: After completing this paper we became aware of similar work on the leading 
N c mass spectrum of excited baryons done in [2^ . 
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APPENDIX A: COEFFICIENTS 

We give here the explicit results for the coefficients Ci_ 8 (9m, @m) used in Sec. III. We 
denote in these expressions t 1 = tan 6^ and t 3 = tan 9^ 3 . 

C 2 (9 N1 ,9 N3 ) = - * - {(1 + tl)(-2tl + t 1 )iV(1650) + (1 + t\)(V^h + 2t§)iV(1700) 

-(1 + t\)(-2 + V^Otg) iV(1520) - (2 + + ^)AT(1535)} (Al) 

C 3 (9 m ,9 m ) = - * - {5(1 + 4tx)(l + tjj)JY(1650) 

+ (-8 - 4 V / T0t 3 )(l + t?)iV(1700) + (4v / 10t 3 - 8t 2 3 )(l + t?)iV(1520) 

+(-20ti + 5t?)(l + ^)JV(1535) + 3(1 + t\)(l + ^)JV(1675)} (A2) 
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Ci(O m ,O m ) = - - - {5(1 - 4$(1 + $#(1650) 

-9(1 + $(1 + $#(1675) + 4(1 + $(1 + 5$#(1700) 

+5(-4 + $(1 + $#(1535) + 4(1 + $(5 + $#(1520)} (A3) 

C 5 (9 m ,9 m ) = - - 1 - {5(3 + 2 tl - 4$ (1 + $#(1650) 

-27(1 + $(1 + $#(1675) + (1 + $(12 + 10VT0t 3 + 20$#(1700) (A4) 
+2(10 - hVWlh + 6$(1 + $#(1520) + 5(1 + *§)(— 4 - 2ti + 3$#(1535)} 

Ct(e m ,0m) + GCiWm, dm) = g(1 - ^ (1 + g {(2 + *?)(! + $#(1535) 
+(1 + 2^)(1 + $#(1650) + 3(1 + $(1 + $#(1675) 

+2(1 + 2$(1 + t?)#(1700) + 2(2 + $(1 + $#(1520)} (A5) 
C 7 (9 m ,e m ) - 6d(g m , 9 m ) = - ^ - {(-4 + $(1 + ^)#(1535) 
+ (1 - 4t?)(l + ^)#(1650) + 3(1 + $(1 + $#(1675) 

+2(1 - 4$(1 + $#(1700) + 2(-4 + $(1 + $#(1520)} (A6) 

c s (e m ,e m ) = - - 1 - {(40 + 8v / Tot 3 )(i + ^)#(i7oo) 

+ (40t! - 25$(1 + $#(1535) + (1 + ^)(-8v / 10t 3 + 40^)#(1520) 

-5(5 + 8ti)(l + tg)#(1650) - 15(1 + $(1 + ^)#(1675)} . (A7) 
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